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Abstract 

The Schrodinger equation for the four-dimensional double singular oscillator is separable in Eule- 
rian, doble polar and spheroidal coordinates in IR 4 . It is shown that the coefficients for the expansion 
of double polar basis in terms of the Eulerian basis can be expressed through the Clebsch-Gordan 
coefficients of the group 5(7(2) analytically continued to real values of their arguments. The coeffi- 

£ — , cients for the expansions of the spheroidal basis in terms of the Eulerian and double polar bases are 

f*"^ ■ proved to satisfy three-term recursion relations. 

o 

1 Introduction 

The Schrodinger equation for a four-dimensional double singular oscillator has the form pQ 

where Ui (i = 0, 1, 2, 3) are Cartesian coordinates of the space IR 4 and c\ and c 2 nonnegative constants. 

In our paper pQ it is shown that the four-dimensional double singular oscillator and generakized 
MIC-Kepler problem decribed by the Schrodinger equation [2] 
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are dual. Here the vector potential A corresponds to the Dirac monopole with the magnetic charge 
g = Tics I e (s = 0, ±1/2, ±1, . . .) and has the form 
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The transformation of duality is the generalization version of the so-called Kustaanheimo-Stiefel trans- 
formation 



CO 

x + iy — 2 (u + iu\) (u 2 + iu 3 ) , 

z = ul + u\ - u\ - u\ (3) 



1 (uo- iui) (u 2 + iu 3 ) 
7 = — In -— - 

2 {uo + iu\) (u 2 - iu3) 



supplemented with the ansatz s — » —id/dj. We also noted, that the parameters of these systems are 
connected with each other by the relations: 

2 

E = 4e 2 , e=—^—, c a = 2A a , where a = 1,2. 
8 

The first two lines of ([3]) are the transformation IR 4 — > H 3 suggested by Kustaanheimo and Stiefel for 
regularization of the equations of the celestial mechanics [3 ■ Later, this transformation found other 
applications as well [4j [5] . The generalized Kustaanheimo-Stiefel transformation ([3]) (the Kustaanheimo- 
Stiefel transformation supplemented with the angle 7) was used for " synthesis" of the charge-dyon system 
from the four-dimensional isotropic oscillator [B]. 

It should be noted that equation @ for A a = and s / reduces tj the Schrodinger equation for 
the MIC-Kepler system [S]. At s = 0, the Schrodinger equation @ is reduced to the Schrodinger 
equation for the generalized Kepler-Coulomb system [9]. In the case when s — and Ai = A2 ^ 0, 
Eq. ^ reduces to the system suggested by Hartmann, at was used for explanation of the spectrum of 
the benzene molecule QUI EJ HU ■ 
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2 Eulerian and Double Polar Bases 



Determine the Eulerian coordinates in 1R as follows: 



(3 i^ti 

uq + iu\ = u cos — e 2 , 





U2 + iuj, = u sin — e l 2 , 



(4) 



where u £ [0, oo), a G [0,27r), /3 G [0,7r], 7 £ [0,47r). In these coordinates the differential elements of 
length and volume, and the Laplace operator have the form 
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if 
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-—- - 2 cos P + —r 



dad'y dj 2 



is the square of the angular momentum operator. 

The Schrodingcr equation fl} in the Eulerian coordinates ((4]) may be solved by seeking a wave function 
ip of the form 

^{u,a,p rf ) = R{u) Zia^P^). (5) 
finding the eigenfunctions of A, J3, J3 j of commuting operators, where the constant 



This amouns to 
of motion A reads 



A = J 2 + — 

h 2 V 1 + cos P 1 - cos ^ 

and which in the Cartesian coordinates Ui has the form 
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The operators J3 and J3 are defined as follows 
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2^ 2 V M o + u i ' u l+ u i 
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(7) 



dtp 
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mip, J 3 ^ 



dtp 
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After substitution the expression ([5]) into the Eq. |T]) the variables in the Schrodinger equation ([T]) 
are separated and we arrive at the following system of coupled differential equations: 



1 d ( . dz_ 

sin P dp \ S1U ^ dp 
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2(1 + cos P) 
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da d"f 
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(8) 
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(9) 
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where A is the separation constant which is the eigenvalue of the operator © . 
The solution of Eq. (jHJ) is easily found to be 

Z jms (a,(3, T ,Si,S 2 ) = N jms (61,62) (cosQ * (sinQ ' P^^\co S p)e ima e is \ (10) 

where mi^ = \m ± s\ + 6i t 2 — *J (m ± s) 2 + rn+ = (\m + s\ + \m — s\)/2 and P r [ a ' b \x) denotes 

a Jacobi polynomial. The quantum number j characterizes the total angular momentum and for the 
(half)integer j the quantum numbers m and s are (half)integer. At a fixed value j the m and s run 
through values: m, s = —j, — j + 1, . . . , j — 1, j. 

Farthemore, the separation constant A is quantized as 

x-( j+ h+>i)( i + h+h + 1 ). („) 



The normalization constant Nj m (6±, 6 2 ) in (JTUJ) is given (up to a phase factor) by 



N jms (6 1 ,S 2 ) = (-l) m " s+ 2 |m " s 



(2j + ^ + 5 2 + 2)(j - m+)\T(j + m+ + 5i + 5 2 + 1) 



y 167r 2 r(j + m_+5i + l)r(j-m_+(52 + l) 
where m_ = (|m + s| — \m — s\)/2 and we assume that 

\JJ J sin P Z j'7n's' {a,f3,j;6 1 ,S 2 ) Z jms (a,P,j]S 1 ,S 2 )dad(3dj = Sjj'5 mm '6 ss/ . 
000 

Let us go now to the radial equation @ . The introduction of (jTTJ) into @ leads to 

^^f)-^W + * + *)W + * + * + ^ + ^(*-^)i»-a 

The solution of this equation normalized by the condition 

u 3 R Nj (u; 61,63) Rn'j (u;6i,6 2 )du = 5 NN , (12) 



has the form 

R Nj (u; 61,62) = C Nj (6i,6 2 ) (au) 2j+Sl+S2 e~^F (~ + j; 2j + 61 + 6 2 + 2; a 2 uA , (13) 



where F(a; c; x) is the confluent hypergeometric function, a = y^Mv/h and 



r (6 6)- Aa " rjJ+j + 61 + 62 + 2) 

C »^M - T(2j + 6 1+ 6 2 + 2^ (fZflj 
The energy spectrum has the form 

E N =hu (N + S 1 +5 2 + 2) (14) 

where N is the principle quantum number, N = 0, 1, 2, . . ., and the quantum number j run through the 
values: j — m + , m + + 1, . . . , N/2. 

Thus, the Eulerian basis (O is the eigenfunction of the operator ([6]) and 

klpNjms = (j + Sl ~t ^ ) (j + Sl ~t ^ + 1 J ^Njms- (15) 
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In the limiting case 8\ = <5 2 = we recover the familiar results for the four-dimensional isotropic 
oscillator [T4"] . 

Let us consider the four-dimensional double singular oscillator in the double polar coordinates [14] . 
In the double polar coordinates pi, p 2 6 [0, 00), <pi,(f2 G [0, 2tt), defined by the formulae 

uq + iu\ = px e llpl , u 2 + 1U3 — p2 e* V2 , (16) 

the differential elements of length and volume read 

dl 2 = dpi + dpi + p\ dip\ + pi dcp 2 , dV = pi p 2 dpi dip 2 , 

while the Laplace operator looks like 

d 2 1 B ( B \ 1 B ( 3 \ Id 2 Id 2 



du 2 pi dpi \ Pl dpi) p 2 dp 2 \dp2J ' p\ dipl ' p\ d(p\' 
The substitution 

where Mi, M2 = 0, ±1, ±2, . . ., separates the variables in the Schrodinger equation and we arrive at 
the following system of equations: 



LA. Ufi) _ ( W W + % ( f + * n - ggT| *, = 0, (i7) 

pi dpi y dpi ) p\ h 2 \ 2 4/x 1 1 



1 d ( d$ 2 \ (|Af 2 | + A 2 ) 2 2p (E h 2 ^p 2 \ 
72diA P2 ^) ?2 + FU"V $2 



(18) 



where f3 - is the separation constant and A a — y M% + 2pc a /h 2 — \M a \. 

These equations are analogous with the equations of the circular oscillator in the polar coordinates 
[To] . Thus, we get 

^n 1 n 2 m 1 m 2 (Pi,P2^i^2;Si,S 2 ) = ^ NlMl (a 2 p 2 ;Ai)^ N2 M 2 (a 2 p2,^2)e lMlV1 e lM ^\ (19) 

where 



, ^ 2r(jV a + |M a | + A a + l) ae-f a;(|M |+A )/2 ■ » L1 , 

= V (Agl r(|M a | + A Q + i) F{ ~ N - |Ma| 1; *>- 

Here N\ and iV 2 are nonnegative integers 

|Afi| + Ai + l E D_ [M 2 | + A 2 + l E _n_ 

1 2 4?iw 8a 2 ' 2 2 4?iw 8a 2 ' 



From the last relations, taking into account (| 14[> . we get that the double polar quantum numbers Ni and 
N2 are connected with the principal quantum number N as follows: 

^ = 2JVi+2JVa + |Mi| + |M 2 |. 
Excluding the energy E from Eqs. (fl7|) and (fl8|) . we obtain the additional integral of motion 



n 



15 / 3 \ I 3 f 3 \ Id 2 13 



P2 dp 2 V 2 dp 2 J pi dpi V 1 9pi y ' (9<p§. p 2 dipl 

(20) 
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with the eigenvalues 



fl=^ (2Nt - 2N 2 - \M X \ + \M 2 \ — Ai + 5 2 ) 
n 

and eigenfunctions ipN 1 N 2 M 1 M 2 (pi, P2, </?i, V2; Ai, A 2 ), i.e. 

Qlp NlN2 Mi M 2 (Pi i />2, ¥>1 1 ¥2\ Ai , A 2 ) = fi^ATj iV 2 Mi m 2 (pi i P2 , ¥> i , <^2 ; A x , A 2 ) . 
In the Cartesian coordinates, the operator Ct can be rewritten as 



d 2 d 2 d 2 d 2 
du\ du 2 8uq du\ 



+ 



+ 



I" u t.,2 , „,2 „,2 „,2\ , 2 ^ ( c l 



(21) 
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3 Connection Between Eulerian and Double Polar Bases 

At fixed energy values we write down the double polar bound states (1191) as a coherent quantum mixture 
of Eulerian bound states 



1pN 1 N 2 M 1 M 2 - J! ^NiN 2 M 1 M 2 



Njms ■ 



(23) 



Our goal is the derivation of an explicit form of the coefficient ^jv™jv 2 MiM 2 • First, we should like to note 
that from the comparison of J4]) with (p~6j) we have 



pi = u cos • 



p 2 = usm-, 



a + 7 



</?2 



a — 7 



(24) 



In relation (|23[) . according to we pass from the double polar coordinates to the Eulerian ones. Then, 
substituting (3 = 0, taking account of 

P M) (1)= (i±^« 

TV. 

and using the orthogonality condition for radial wave functions in hypermomentum [14j 
Rjqy (u; S%, S 2 ) Rnj (u; Si, S 2 ) udu 



2a 2 5 jr 



(2j + St + S 2 + 2) ' 



we obtain the following integral representation for the coefficients M 



M 2 



W, 



pns 
NiN; 



y/{2j + S 1 +S 2 + l){j-m+)\ 



Mi m 2 - l^ 2+ X)T(2^ + S 1 +S 2 '+2) E °™ N > K "™ SMl ' m+s SaL - 



Here 



Cj N l N 2 



; r [ ^+j + s 1 + s 2 + 2) ■ 



r (j - m_ + S x + 1) r(A^! + mi + l)r(iV 2 + m 2 + 1) 
(iV 1 )!(AT 2 )!(f - j)W (j + m_ + S 2 + 1) T (j + m+ + S x + S 2 + 1) 



(25) 



(26) 
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and 



K 



Nm 



e ~x x j+m 1 +m 2 +8 1 +5.; 



Fl-Xr.un + l:.r) F [ -^+j;2j + S 1 +S 2 + 2;x) ,1,-. (27 1 



where x = a 2 u 2 . Further, in (|27p writing down the confluent hypergeometric function F (— N±; mi + 1; x) 
as a polynomial, integrating by the formula [16] 

oo 

J e~ Xx x u F (a; 7; kx) dx = ^-^5^2-^1 (a: v + l',7^ dx, 



and taking account of the relation 



2 F 1 (a, b; c; 1) 



T(c)T(c-a-b) 
T(c-a)T{c-b) 



we derive 



NNl = (f - m + )\T(2j + Si + S 2 + 2)r (j + m + + ft + 6 2 + 1) 
(i-m + )!r(f +j + S 1 + S 2 + l) 



x 3 F 



-Nx, -j + m + ,j + m+ + Si + S 2 + 1 
2<] mi + l,-f + 1 



The introduction of ((26j) and (J25J) into ((25j) gives 



(2j + gi + ft + 1) r(jVx + mi + l)r(JV 2 + m 2 + 1) 
(JVi)!(JV 2 )! (f - - m+y.T (j + m_ + ft + 1) 



(t -"■»+)! / r (j - m_ + ft + 1) r {j + m + + ft + ft + 1) 
r(mi + l) V 



r(f + j + ft + ft + i) 
1} tfj* 

It is known that the Clebsch-Gordan coefficients for the group SU (2) can be written as [9] 

(2c + 1)(6 - a + c)l(a + a)\(b + (3)\{c + 7)! 



-m,-j + m + ,j + TO + + ft + ft + 1 
X3F2 ^ m 1 + l,-f +m + + l 



,--rC7 



(b - /3)!(c - 7)!(o + & - c)!(a - b + c)!(a + 6 + c + 1)! 



1/2 



vA (~l) a a (a + 6-7)! f -a + a,c + 7+l,-c + 7 
7 ' a+/3 V^^)!(6-a + 7)! 3 2 I 7 -a-&,&-a + 7 + l 

Finally, comparing (|2T)|) and (f5TT|) . we arrive at the following representation: 



Ni- 



SMi,m+s SM 2 ,m—s x 



xC" 



f + 2m_+25 2 -2 7ji 2 + Af 2 -7V 1 JV -2m _ +25-1 - 2 r^H-Ni— _Wo 



(28) 



(29) 



(30) 



(31) 



Equation (|3"Tj) proves that the coefficients for the expansion of the parabolic basis in terms of the spherical 
basis are nothing but the analytical continuation, for real values of their arguments, of the SU(2) Clebsch- 
Gordan coefficients. 
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The inverse representation has the form 



4>Njms — ^ W^ ms (5ii6 2 ) 1pN 1 N 2 M 1 M 2 - 
Ni,MiM 2 



(32) 



The expansion coefficients in (|32l) are given by the expression 



Mi,m+s 0M 2 ,m-s 



xC 



3 + 



*l+ 5 2 rn 1 +m 2 



N + 2 tti_+ 2 S 2 — 2 N + 2 m_+ 2 6 2 - 2 N~ 2 jn_+ 2 S 1 — 2 . N- 2 m_~ 2 S 1 - 2 - 

1 , j "i; 4 ,ni+\m-s\ 4 



(33) 
(34) 



4 Prolate Spheroidal Basis 

Let us determine the four-dimensional spheroidal coordinates 



u a +iu x = - v / (C + l)(l + ^)' 



u 2 + iu=- \/(£, - - v) ' 



(35) 



where £ 6 [1, oo), r\ G [— 1 , f], and d is the interfocus distance. 

In the spheroidal system of coordinates the four-dimensional double singular oscillator potential has 
the form 



C-2 



In the coordinates (|35|) Laplace operator have the form 



( 9 2 



8 [9 



(e - 1) h 



1 / 1 



I 



d d 



1/1 



9 " 






dr] 
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1 



,9 a 



2VC + 1 1 + 17/ \da djj 2 V£ - 1 l-77/V<9a<97 
After the substitution 

the variables in the Schrodinger equation |T]) are separated 



d „\ e£ m? m 2 , 

— (P - 1) 1 - — 

d£^ >dC 2(^ + 1) 2(C - 1) 



16 

a 4 d 4 



— (l - 77 2 ) 

drj K 11 dr] 2(1 + 77) 2(l-r?) 16 



(e - 1) + ^ 



2 , fj,Ed 2 1 



'02 = —Qlp2, 



(36) 
(37) 



where Q is a separation constant in spheroidal coordinates. By eliminating the the energy E from Eqs. (|36 
and (|37p . we produce the operator 



n 1 I 9 

£-77 [ <9£ 

e + 1? + 1 



d d 



dr] 



( ^ d-q 



£ + r?-l 



+ 



2 (£ + 1) (1 + 77) V<9a ' d~fj 2 (£ - 1) (1 - rj) \da ~ !fy 
ci^-i?) , c 2 {^-ri) 



(38) 



4 j4 

a a 



(^ + 1) 



2(£+l)(l + i?) 2(C-l)(l-r / )' 
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the eigenvalues of which are Q and the eigenfunctions of which are ?/;(£, 77, a, 7). The significance of the 
self-adjoint operator Q can be found by switching to Cartesian coordinates. Passing to the Cartesian 
coordinates in ([55)) and taking (0 and (J^U) into account, we obtain 



Therefore, 



a 2 d 2 - 
4 



QlpNqms(£, V, Q> 75 d, £l> fo) = Qq^Nqmsi^, V, a, 7; rf, (Si , <5 2 ) , 



(39) 



(40) 



where index q labels the eigenvalues of the operator Q. 

Now construct the spheroidal basis of the four-dimensional double singular oscillator using the follow- 
ing expressions: 



N/2 



^Nqms = U NqmA R i 5 ^ S ^) ^ 



Njms • 



(41) 



(N-\M 1 \-\M 2 \)/2 



(42) 



JVi=0 

Substituting (|4Tj) and (|42|) into (|40)l . and then using (|39|) we arrive at the following algebraic equations: 

a 2 d 2 



TP _ 

^ Nqms 1 



a 2 <i 2 / mi — mo 

ATi - N 2 + 



A'; 



where 

(Q) jr = J 1p* Njms Cll/j Njms dV, = J TpN 1 N 2 M 1 M 2 ^i } N[N! 2 M 1 M 2 dV. 

Now using expansions ([3"T|) . (|33|l and formulae [17] 

4c 2 (2c+ l)(2c- 1) 



(~IC1 _ 



-,1/2 



(c + 7) (c - 7) (b - a + c) (a - b + c) {a + b - c + 1) (a + b + c + 1) _ 
(c - 7 - l)(c + 7 - 1)(6 -a + c- l)(a -b + c- l)(a + 6 - c + 2) (a + 6 + c)l 1/2 



4(c- l) 2 (2c-3)(2c- 1) 



vriC -2. 7 _ (a - /?)c(c - 1) - 70(0 + 1) + 76(6 + 1) ^ c _i, 7 I 

XO aa;fc/3 2c(c - 1) V¥ [' 



(43) 



[ c(c + 1) - a(a + 1) - 6(6 + 1) - 2a/3] C a c : 2; ;fci/3 = 

= VO* + «) (a - « + - + /3 + 1) <^_ 1;W + 
+ v /(a-a)(a + a + l)(b + /3)(b-0 + l) C c a l +Vb ^_ v 



and with the ortho-normalization conditions 



a+b 



7'7> 



E^C7 sycy _ r 



c— 1 7 



for the Clebsch-Gordan coefficients of the group SU(2), for the matrix elements (S2) 3 y and (A^at; wc 
get the expressions 



,Ax . = (mi +m 2 )(mi -m 2 ) _ 2 (A^fyj+i + A^fyj.i) 



(2j + <5i + «5 2 )(2j + tfi + <5 2 + 2) 



2j + + <5 2 



(44) 



where 



(A) Wl Ar; 



iV 

(JVi + 1)(AT 2 + m_) + (— - ATi + <5 2 )(iVi + |m - s| + fe) + 



+m_(m+ + 5 2 ) -(^i - o" 2 )(o"i - £ 2 - 2)+ 



Ai 



Nms 



~VMNi + l)(JVi + mi + l)(iV 2 + m 2 )6 N ^ Nl+1 - 
-y/N^Ni + l)(JVi + mi + 1)(JV 2 + m 2 + 1)^,^. 



V (j - m_ + <5i)(j + m_ + <5 2 ) x 

C? - + m + + Si + ggKjV - 2j)(N + 2j + 2fr + 25 2 ) 

4 ^ + 5i±fe) 2 (2j - + ^ + j 2 _ i)( 2 j + s t + S 2 + 1) 



1/2 



(45) 



Substituting expressions (|44|) and (J45f into the algebraic equations (|43|) , we derive the three-term recursion 
relations 



4J+1 7-/3 + 1 4j XjJ-l 

-^Nms Nqms ~ ^-nms u Nnqms 



2 d 2 



J H ^ — L? + — z 1- 1 



(mi + m 2 )(mi - m 2 ) [ ^ 



(2j + <5i + <5 2 ) (2j + 5i + <5 2 + 2) [ N « ms 



(JVi + l)(JV a + m_) + (JV - Nt + S 2 ){N 1 + m 2 ) + -{8 X - 5 2 )(S 1 - ff 2 - 2) + 



. . . a d / r „ T mi — m 2 . 

+m_(m+ + <5 2 ) + — f Ni - N 2 + - ) - Q q 



V, 



Ni 



Nqms 



= VMNi + lWi + ™i + l)(N 2 + m 2 ) V^ q +) + 
+ VNi(N 2 + l)(JVi + mi + 1)(JV 2 + m 2 + 1) V$) ~] 
for the interbasis expansions coefficients U 3 Nqms and V^ q 1 ms . 
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